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By introducing the thermo entangled state representation, we convert the calculation of Wigner 
function (WF) of density operator to an overlap between "two pure" states in a two-mode enlarged 
Fock space. Furthermore, we derive a new WF evolution formula of any initial state in self-Kerr 
Medium with photon loss and find that the photon number distribution for any initial state is 
independent of the coupling factor with Kerr Medium, where the number state is not affected by 
the Kerr nonlinearity and evolves into a density operator of binomial distribution. 
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I. INTRODUCTION 

Nonclassicality of optical fields has been a topic of 
great interest in quantum optics and quantum informa- 
tion processing [l[ , which is usually associated with quan- 
tum interference and entanglement. The phase space 
Wigner function (WF) Q, Q of quantum states of light 
is a powerful tool for investigating such nonclassical ef- 
fects. The WF was first introduced by Wigner in 1932 to 
calculate quantum corrections to a classical distribution 
function of a quantum-mechanical system. The partial 
negativity of the WF is indeed a good indication of the 
highly nonclassical character of the state [4] and moni- 
tors a decoherence process of a quantum state, e. g. the 
excited coherent state in both photon-loss and thermal 
channels 0, Q, the single-photon subtracted squeezed 
vacuum state in both amplitude decay and phase damp- 
ing channels 0, and so on (8l-[l2j. 

Nonlinear interaction of light in a medium provides a 
very useful framework to study various nonclassical prop- 
erties of quantum states of radiation. The Kerr medium 
is one of the simplest nonlinearity, which shall allow us 
to investigate the full time-dependent WF dynamics with 
or without a quantum noise. Recently, a Fokker-Planck 
equation for the WF evolution in a noisy Kerr medium 
(^( 3 ) nonlinearity) is presented [l3j|. Then the authors 
numerically solved this equation assuming coherent state 
as an initial condition and discussed its dissipation ef- 
fects. However, for any initial condition, as far as we 
are concerned, there is no report about the WF evolu- 
tion. On the other hand, by using the thermo entan- 
gled state representation (TESR) we solved various mas- 
ter equations to obtain density operators with an infinite 
operator-sum representation [14[ and then revealed that 
the WF of density operator can be expressed as an over- 
lap between two pure states (see Eq. ([20]) below) [Hj]. 
This brings much convenience to calculate time evolu- 
tion of WFs when quantum decoherence happens. Thus 
the TESR is beneficial to quantum decoherence theory. 

In this paper, we shall appeal the TESR \rj) to treat 
the WF evolution at any initial condition in self-Kerr 



Medium with photon loss and present a new formula to 
calculate time evolution of the WF for quantum deco- 
herence. In addition, based on the derived WF evolution 
formula, we shall deduce the photon number distribu- 
tion for any initial state in presence of Kerr interaction, 
where the photon number distribution is independent of 
the coupling factor \ that is relative to the Kerr medium. 
As examples, the WF formula is applied to the cases 
of initial coherent state, and number state, respectively. 
Conclusions are involved in the last section. 



II. BRIEF REVIEW OF THERMO ENTANGLED 
STATE REPRESENTATION 

We begin with briefly reviewing the thermo entangled 
state representation (TESR). On the basis of Umezawa- 
Takahash thermo field dynamics (TFD) [l6l4l8[ we con- 
structed the TESR in the doubled Fock space[l9|-[2l| , 
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is the displacement operator. 
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is a fictitious mode accompanying the real photon cre- 
ation operator aft , |0,0) = |0) |5) , and |5) is annihilated 
by a with the relations [a, aft] = 1 and [a, aft] = 0. The 
structure of \rj} is similar to that of the EPR eigenstate 
shown in Ref. [l9|. Operating a and a on \r]) in Eq.(pQ) 
we can obtain the eigen-equations of |^), 
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Note that [(a — eft), (aft — a)] = 0, thus I77) is the 
common eigenvector of (a — eft) and (a — aft). Us- 
ing the normally ordered form of vacuum projector 
|0,5) (0,5| =: exp (—aft a — aft a) : and the technique of 
integration within an ordered product (IWOP) of opera- 
tors [221-1241]. we can easily prove that I77) is complete and 
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orthonormal, 
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It is easily seen that \rj = 0) has the properties 

CO 
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(where n = n, and n denotes the number in the fictitious 
Hilbert space) and 

a \rj = 0) = a ] \rj = 0) , 
a t |r/ = 0) = a I77 = 0) , (5) 
(a t a) n |r/ = 0) = (a+a)" |r/ = 0) . 

Note that density operators p(a^a) are defined in the 
real space which are commutative with operators {a\a) 
in the tilde space. 

In a similar way, we can introduce the state vector |£) 
conjugated to |^), defined as 

|£)=D(£)e- atat |0,0) 

= exp ^l^ + ^+fat-a^j |0,6) 

= (-l)« t «|r ? =-0, (6) 

which also possesses orthonormal and complete proper- 
ties 
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III. MASTER EQUATION FOR A SELF-KERR 
INTERACTION 

In the Markov approximation and interaction picture 
the master equation for a dissipative cavity with Kerr 
medium has the form [25|, [26| 



dp 
~dt 
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where 7 is decaying parameter of the dissipative cavity, 
X is coupling factor depending on the Kerr medium. Mil- 
burn and Holmes [27j solved this equation by changing 
it to a partial differential equation for the Q-function 
and for an initial coherent state. Here we will solve the 
master equation by virtue of the entangled state repre- 
sentation and present the infinite sum representation of 
density operator. 

Operating the both sides of Eq.flSJ) on the state \rj = 0) , 
letting \p) = p\r] = 0) (Here one should understand the 
single- mode density operator p in the left of Eq.® as the 
direct product p<&I when p acts onto the two- mode state 



1 77 = 0) = e atat |0, 0), where / is the identity operator in 
the auxiliary mode.), and using Eq.(|5j) we can convert 
the master equation in Eq. J5J) into the following form, 
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dt 



\p) = j-ix (a j a) 2 - (d ] dy 

+ 7 (2ad — a) a — a) a) } \p) , (9) 



i.e., an evolution equation of state vector \p). Its solution 
is then of the form 



\p) 
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where \po) = po \tj = 0) , po is an initial density operator. 
The advantage of using thermo field notation over more 
traditional algebraic manipulation with superoperators is 
that in many situations (and, particularly, ones of our in- 
terest) it enables to simplify, make more illustrative and 
less cumbersome finding the solution ([TP]) and estima- 
tion of time-dependent matrix elements. In particular, 
it allows to represent in a simple form a factorization of 
the superoperator exp{- • • } into multipliers with easily 
estimated actions on the number states [28j. 
By introducing the following operators, 

+ +~ a ^ a + + 1 
ivo = a 1 a — a 1 a, K z = , A_ = aa, (11) 

which satisfy [Ko,K z ] = [Ko,K-] = 0, we can rewrite 
Eq.$T0j) as 

\ p \ = e {-ixt[Ko{2K z -l)] +1 t{2K_-2K z + l)} |^ 



= exp [ixtK + 7*] 
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With the aid of the operator identity |30| 

e X(A+*B) = e XA exp r aB ^ _ g -Ar) 

= exp [oB (e AT - l) /r] e XA , 



(13) 



which is valid for [A, B] = tB, and noticing [K Z1 K_] = 
— we can reform Eq. ([T2]) as 

\p) = exp [ixtKo + Jt] exp [T Z K Z ] exp [r_if_] |p ) , 

(14) 

where 
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From Eq. ([T4]) we can obtain the infinite operator-sum 
form of p(t), (see Appendix A) 
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where the two operators M m ^ nj i and M m n l are respec- 
tively defined as 




Although M m ^ n j is not hermite conjugate to n j, the 

normalization still holds, (Em, n ,z=o ^ln,n,i M m,n,l = 1, 
see Appendix B |29|) i.e., they are trace-preserving in a 
general sense, so M m?n5 / and n l may be named the 
generalized Kraus operators. 



IV. EVOLUTION OF WIGNER FUNCTION 
FOR SELF-KERR CHANNEL 

In this section, we consider Wigner function's time evo- 
lution in the self-Kerr medium channel. For this purpose, 
we shall derive a new expression of Wigner function in the 
TESR. According to the definition of Wigner function of 
density operator p, 



W(a,a*) = Tr [A 



(18) 



where A (a, a*) is the single- mode Wigner operator |2|, 
|30| , whose explicit normally ordered form is (3lj 



A (a, a*) 



7T 7T 

(19) 

By using (h\ in) = 5 ny7n (n = n, m = fa) and noticing (j6j) 
as well as \p) = p \rj = 0) , we can reform Eq. ([T8]) as 
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where Eq. (|20j) is the Wigner function formula in thermo 
entangled state representation, with which the Wigner 
function of density operator is simplified as an overlap 
between two "pure states" in enlarged Fock space, rather 
than using ensemble average in the system- mode space. 
This will brings much convenience to calculate the time 
evolution of Wigner functions when quantum decoher- 
ence happens. 

Projecting (fTl|) on the entangled state representation 
- (£=2a| ? an d inserting the completeness relation flTJ), we 
find 



W(a,a*,i) =4 



J— G(a, 



/M)W(/3,/T,0), (21) 



where W (a, a*,t) and W (/?, /?*, 0) are the Wigner func- 
tions at the evolving time t and initial time, respectively, 
and 

G(a,0,t) = (£=2«\exp[ixtK + 'Yt] 

x exp [Y Z K Z ] exp [r_K_] |C=2/3> • (22) 



It is convenient to the matrix element in ([22]) according 
to the two-mode Fock space. Thus the (£=2a| is expanded 
as 
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By using the two-mode Fock state |m,n) 
a |m fl tn /vS! |0, 0), we get 



(e|m,n)=if m , n (r,0e- |el /2 /V^!, (24) 

where iJ m.n ( £* 2 £ ) is the two- variable Hermite poly- 
nomials [32|, |33j. Inserting the complete relation 
n=o l m > ™) ( m > ^1 = 1j after a long but straight calcu- 
lation, then the Wigner function's evolution is given by 
(see Appendix C) 
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It is obvious that, when x — 0, the case of photon loss, 
A m , n (1 - e~ 2lt ) = T and Eq.(j25]) just does reduce to 
(see Appendix E) 

W(a,a*,t) = ^J ^c-*M^lV03,/r > 0) > 

(28) 

which is just the evolving formula of Wigner function for 
amplitude-damping channel. While for 7 = 0, without 
photon- loss, Eq. (|25]) reduces to 

W(a,a*,i) 

_ ~ exp [-itf (m 2 -n 2 )] 
~ 1^ , , oui2 H my7l (2a ,2a) 
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V. PHOTON NUMBER DISTRIBUTION IN 
PRESENCE OF KERR INTERACTION 

Now we consider photon number (PN) distribution in 
presence of Kerr medium. According to the TFD, we can 
reform the PN p (n) = tr [p \n) (n\] as 

oo 

p(ri) = (n\p\n) = (n,n\p\m,fh) 



m=0 



(n,h\p\r] = 0) = (n,n| p) , 



(30) 



which is converted to the matrix element (n, h\ p) in the 
context of t her mo dynamics. Then using the complete- 
ness of (£| and Eq. ([3Q|) as well as Eq. (|2Q|) . we have 

P(n) = J ^<n,n| f) (£| p) 

= J d 2 £ <n, h\ £)W(a = e/2, a* = £*/2) 

= 4^ d 2 aW ln){nl (a, a*) W (a, a*) , (31) 

where (a,a*) = ^e' 2 ^ L n (4 \a\ 2 ) is the 

Wig ner function of number state \n) (n| as shown in 
(34 [35j. Thus one can calculate the PN by combining 
Eqs.(|2S]) and (|3T|) . 

Next we evaluate the PN for the above decoherence 
model in Eq.flHJ). Substituting Eq.(|25]) into Eq.flU]), we 
have 

00 ^ 
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Then substituting Eqs.(j33]) and (|27j) into (|32j) yields 
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VI. EVOLUTION OF QUANTUM STATES 

The phase space Wigner distribution function descrip- 
tion of quantum states of light is a powerful tool to inves- 
tigate nonclassical effects, such as quantum interference 
and entanglement. In this section, as the applications of 
WF evolution formula, we take two special initial states 
as examples. 

(1) When the initial state is the coherent state \z), 
whose WF is given W (/3,/3*) ± e ~ 2|/3_2;| , thus substi- 
tuting into Eq. (|27|) yields (see Appendix G) 
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which is a new expression of the evolution of WF for 
any initial state. In particular, when 7 = 0, without the 
dissipation, Eq. (|36|) reduces to 



- n2 )ff m>n (2a*,2c$,7) 
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which is identical to Eq.(7) in Ref.[l3j], where Eq.(7) is 
used to make the numerical calculation since it is much 
more rapid than the other expression Eq.(6). In addition, 
from Eq. (|36|) one can see that the WF can be obtained 
very quickly when the dissipation cannot be negligible. 
Further when xt = 27r, Eq. ([37|) just returns to the WF 
of the initial coherent state. 

Fig.l presents the plots of the WF for different pa- 
rameters xt an d a = 2. From Fig.l, one can see that 
the WF turns into an ellipse and squeezing appears in 
an appropriate direction. Then the ellipse changes into 
a banana shape and a tailor of the interference fringes 
appears where the distribution takes the negative values. 

(2) Another example is number state, where the WF 
of number state \s) is given by 



W a (P,l3*,0) 



1/sl 

7T 

(-1)* 



' 2 ^ 2 H StS (2/3,2/3*) 
e- 2 ^l 2 L s (4|/3| 2 ), (38) 



which corresponds to the photon number of density op- 
erator in the amplitude-damping quantum channel [15| . 
From Eq. ([34]h it is easily to see that, for any initial state, 
the photon number distribution p (s) is independent of 
the coupling factor x that is relative to the Kerr medium, 
as respected in [36|. 



substituting it into (|27|) and using the generating function 
of two- variable Hermite polynomials [see below Eq.(F3)], 
we have 
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FIG. 1: (Color online) The WF for different parameters 
and a = 2, (a) # 0, (b) x t 0.04, (c) x* 0.06, (d) 
0.08, (e) x* = 0.1, (f ) x* = 0.2. 



thus the evolution of WF for \s) is 

W a (a,a*,t) 

-2<yt\ s - rn 



ra=0 x 



W m (a, a*, 0X40) 



which is the WF of number state in the photon-loss chan- 
nel and indicates that the number state is not affected 
by the Kerr nonlinear ity. In particular, when 7 = 0, or 
t = 0, Eg. (pTO]) just reduces to the WF of number state. 
From Eq. ([40|h on the other hand, it is found that the 
number state evolves into a density operator of binomial 
distribution (a mixed state) if e~ 2lt is a binomial param- 
eter. 



VII. CONCLUSIONS 

In summary, by converting Wigner function for quan- 
tum state into an overlap between two "pure states" in 
a two- mode enlarged Fock space, we investigate the WF 
evolution of any initial condition in self-Kerr Medium 
with photon loss and present a new formula for calculat- 
ing time evolution of the WF for quantum decoherence. 
Based on the derived WF evolution formula, in addition, 
we discuss the photon number distribution for any initial 
state in presence of Kerr interaction. It is found that the 



photon number distribution is independent of the cou- 
pling factor x i n correlation with the Kerr medium, as 
expected by people. As applications, furthermore, the 
two cases of initial coherent state and number state are 
considered. It is shown that the coherent state can be 
squeezed due to the presence of Kerr medium, while the 
number state is not affected by the Kerr nonlinearity and 
evolves into a density operator of binomial distribution 
(a mixed state) with e~ 2lt being a binomial parameter. 
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Appendix A: Derivation of Eq. ([16p 

In order to obtain the infinite operator-sum form of 
p(t) from Eq. (fT4]h using the completeness relation of 
Fock state in the enlarged space ^ n=0 |m,n) (ra, n\ = 

1 and noticing aJ l \n) = \J \n-\-l), we have 



| p ) = e ixt*o+7t e r.*. e r_K_ |po) 
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Furthermore, using the relations 

(n\ 77 = 0) = \n) , |ra, n) = \m) (n\ 77 = 0) , 
we find 
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(ra, n\ a l p a^ 1 \rj = 0) = (m\ (h\ a 1 p Q o) 1 \rj = 0) 
= (m\a l p ^ l ((n\rj = 0)) 
= (m\ a 1 pqo) 1 \n) . (A4) 

Thus Eq.(Al) becomes 
A' 



|^ _ m ' n c -ixt(m 2 -n 2 )-ft(m+n) 

m. ,n. /=0 

x \m, h) (m\ a l poa) 1 \n) 
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x e - ix t(m*-n*)-^ m+ n) \ m ^) p 0Mn+l , (A5) 
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where po,m+z,n+z = (m + l\ po\n + I) . Using Eq.(A3) After depriving \rj = 0) from the both sides of Eq.(A6), 
again, we see the solution of master equation (|8]) appears as an infinite 

operator-sum form 



~ V(n + Q!(m + Q! 

x e -i xf ( m ^)- 7f(m+ „)^ ro+; ^ lm) (nl ^ = 0) _ 

(A6) 



- / („ + Ql (m + Ql A^ e - ixtK _^)- 7t(TO+7t) + + 

^— ' V n!m! /! 



m,n,l=0 



(A7) 
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Note that the factor (m — n) appears in the denominator 
of A m)U (see Eq.(A2)), (this is originated from the non- 
linear term of (o)a) 2 ) so that it is impossible to move all 
n— dependent terms to the right of a 1 p^a) 1 . Fortunately, 
we can formally express Eq.(A7) as Eq. ([T6]) . 

Appendix B: Proof of normalization for the gen- 
eralized Kraus operators 

Using the operator identity 

e Aata =: exp [(e A - l) a+a] : and the IWOP tech- 
nique, we can prove that 



from which one can see that the normalization still holds, 
i.e., they are trace-preserving in a general sense, so 
M m>nj ; and n l may be named the generalized Kraus 
operators. 
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n) (n\ a 1 
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-2i 7 y 



■: exp [(e" 27 * - l) a^a] (a^a) 



= 1, 
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Appendix C: Derivation of Eq. ([25|) 



Using Eq.flU]) and (Al) as well as (A2), Eq.(j22j) can 
be rewritten as 



G(a,(3,t) = e -^(™ 2 -» 2 )-7t( m +n) ^ ^ (m> ft | e A m , re aa |^ 



-2|a| 2 00 A ^ 
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Further using a new sum formula (see appendix D) where C m , n (a,a*,t) is defined in ([26]) . Substituting 

Eq.(C3) into flZQ) yields (p5j) and (p7j). 



^2—H m+ i tn+ i (x,y) 
1=0 L 



Appendix D: Derivation of Eq.(C2) 



Hrn,n \ , . , > = . ) , (C2) 



( z _|_ -g( m + n + 2 )/ 2 m ' n \ + 1 ' + 1 y Using the integratal expression of two-mode Hermite 

polynomials, 

thus Eq.(Cl) can be recast into the following form 



G(a,/3,t)= ^2 C m ,n(a,a*,t)e 



m,n—0 v. f Z 1 1 2 1 /- * 



H m , n (£,r 1 ) = (-l) n e^ f 

x g m ,„ ( -7=^==, /A 2/? * ; ■ ) , (C3) ' " (Dl) 

we have 
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- xz — yz 



\z\ -\- xz — yz 



e «+i 



(a + 1) |z| + xz — yz* 



d 2 z 



(a + l) (m+n+2)/2 



(a + l) (m+n+2)/2 Vv^ + T'v^ + T 



7T 
2/ 
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(D2) 



thus we have completed the proof of (C4). 
Appendix E: Derivation of Eq. ([28|) 

When x = 0, A m , n (1 - e~ 2 ^) = T, and 



Cm,n (ot,Ot*,t) 

_ exp [—jt (m + n)] 
= m!n! (T + i)( m+n + 2 )/ 2 

then we have 



tf m , n (2a*,2a)e- 2 H 2 , (El) 



£ C m , n (a,a*,t)i7 m , n (2/3,2/3*) 



m,n=0 



-2H 2 00 



(x/T+l) 



m+n 



m!n! 



m,n=0 

xff m ,„ (2a*,2a)ff m ,„ (2/3,2/3*)- 



(E2) 



Using the following formula 



—r-rH min (x,y)H min (a,P) 
^-^ mini 

m,n=0 



1 



1 - St 



exp 



sxa + — (xy + a/3) s£ 



Eq.(E2) can be rewritten as 



(E2) 



2T 



Thus Eq. ([25|) becomes 



(E3) 



(E4) 



8 



/,2 o 00 
l{ m,n=0 



2(3 2(3* 



= R.H.S. of Eq.©. 



T- 1 
J T + 1 



.Vr+T' VtTT 

2e~^ (a* (3 + a/3*) - 2e- 2 ^ £ f |a| 2 + 



(M 2 + ^) 



W (/?,/?*, 0) 



(E5) 



Appendix F: Derivation of Eq. ([33|) thus Eq.(F2) becomes the right hand side of Eq. ([33j) . 

Using the relation between Hermite polynomial and 
Lagurre polynomial, 



L m (xy) 



(-1Y 



ml 



(Fl) 



we can recast the left of Eq. (j33p into the following form 
F m,n = ^ J ^e-^ 2 (-l) s s\L s (4\a\ 2 )H m , n (2a*,2a) 
= -. f —e- A ^ 2 H atS (2a*,2a)H m , n (2a*,2a) 

S\ J 7T 

= l/4 r ^ e -\ a f H ■ ( a * a ) H m>n (a*,a). 

S\ J 7T 



(F2) 



Appendix G: Derivation of Eq. (|35p 

For this purpose, from Eq. ([27|) we have 



Further using the generating function of H m n (e, e), 

Qm-\-n 



and the integration formula, 
72 



exp \-tt' + et + | t=t/=0 , (F3) 



^m,n = 4 y ^e- 2 ^^l 2 ^ m , n (2x0,2x0*) e 



■2|/3-*| 2 



/ <L± e W+t z+r > z * = -i e -^,Re(C) < 0, (F4) 



(Gl) 



where we have set 



we have 

q2 



Qm-\-n 0m-\-n 

dt m dt /n dr m dr /n 

i2 



1 A m n 

y = t—x — ^ 



exp \-U' - tt) 



1 + A m>n ' ^ + 1 ' 



(G2) 



exp 



|ar + (t + r) + + r')a 



t=t'=T=T'=0 



d 



m +n Qm-\-n 



dt m 'dt ,n ' dr m dr ,n 
m!nW m / >n (J n / )m , 



exp [tr / + rt / ] £=£/=r=r/=0 



Using Eqs.(F3) and (F4), Eq.(Gl) can be recast into the 
(F5) following form, 



9 



E _ 4p -2M 2 _^ -tt' 



-2M 2 



dt m dt'> 

Qm-\-n 

dt m dt ,n ( 



d 2 f3 



< i - 

J TT 



exp 
T exp 



2 (2/ + 1) |/3| 2 + 2/3 (a* + **) + 2/3* (art' + s) 
|/3| 2 + /3(^ + £*) + /3* (zt' + z) 



t=t'=o 



V + 1 ,0,2 



t=t'=0 



7T i/ + 1 



-e 2/+ 



T M 2 



dt m dt r 



exp 



2x 2 
2/+1 

I — 



2xz 2xz* 




2/+1 J/+1 



t=t'=0 



Note that 1 - ^ = 0,2/ + 1 = 



l-Arr 



3/+1 



2 2?/ — A™ „ - 1, then we find 



l+An 



+ 1 = 



12 __2ilu|2 / 2xz 

/7 1 — y+l 1 1 



2xz* 



7r 2/ + 1 \y + 1/ V^ + :l -, 

- (A m ,„ + l) - ^ e^™-- 1 ^! z m z* n . (G4) 



(G3) 
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